Optics of Nonuniformly Moving Media 
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A moving dielectric appears to light as an effective gravitational field. At low flow velocities the 
dielectric acts on light in the same way as a magnetic field acts on a charged matter wave. We 
develop in detail the geometrical optics of moving dispersionless media. We derive a Hamiltonian 
and a Lagrangian to describe ray propagation. We elucidate how the gravitational and the magnetic 
model of light propagation are related to each other. Finally, we study light propagation around a 
vortex flow. The vortex shows an optical Aharonov-Bohm effect at large distances from the core, 
and, at shorter ranges, the vortex may resemble an optical black hole. 
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I. INTRODUCTION 

Consider a glass container filled with a transparent liq- 
uid, say water. Let a plane wave of coherent laser light 
travel through the water. Obviously, not much will hap- 
pen. The light will remain a plane wave and will only 
gather an overall phase shift. Now imagine that the wa- 
ter is set into motion. For example, a magnetic mixer 
at the bottom of the container creates a vortex. Let us 
assume that no air bubbles contaminate the transparent 
liquid and that no heat gradient is generated. We know 
that water is to a large degree incompressible. Therefore 
the refraction index of the whirling liquid is spatially uni- 
form. Will the light remain a plane wave? Maybe sur- 
prisingly, it will not. Instead, the light will develop an 
interference structure that is sensitive to the velocity of 
the liquid. Furthermore, if we send in a narrow laser 
beam, the vortex will bend the ray. A moving medium 
drags light. This effect can be employed to gather infor- 
mation about the flow of a transparent liquid. One could 
think of reconstructing an unknown velocity profile from 
measured interference patterns, as a form of optical to- 
mography. 

In this paper we develop a systematic theory that 
explains motional effects of a nondispersive dielectric 
medium on light propagation. We postulate that the 
wave equation is valid in all locally comoving frames of 
the medium. Then we transform the wave equation to 
the laboratory frame. In the limit of geometrical optics 
we find the Hamiltonian that determines the trajectories 
of light rays. Using a different implementation of the 
same idea (transforming Maxwell's equations from the 
comoving to the laboratory frame) Berry and Klein |^] 
have also derived effective scalar and vector potentials 
and a Hamiltonian governing light rays and waves. In 
addition, we develop a covariant theory of light propaga- 
tion in moving media. In accordance with earlier papers 
by Gordon |a] and Pham Mau Quan 0] we find that 
light rays follow zero-geodesic lines measured with re- 
spect to a certain curved metric in space-time, similar 
to the light propagation in general relativity. The metric 



of the "glass of water" depends on the refractive index 
and on the fiow, and it establishes a fascinating analogy 
between moving media and gravitational fields. 

Can we see some of the spectacular effects of general 
relativity in an earthly laboratory? Most probably, this 
would take flow velocities that are comparable with the 
speed of light in the medium. Recently, dielectrics with 
incredibly low group velocities have been created [El. 
These media are far from being as simple as ordinary 
liquids, and they are highly dispersive. For instance, the 
refractive index reaches unity, i.e. the phase velocity ap- 
proaches c, at the frequency where the group velocity 
is lowest, i.e. where the refractive index changes most 
rapidly. However, as we show in a separate paper [^, 
many phenomena that are conceivable for dispersionless 
media find an experimentally feasible analog in disper- 
sive dielectrics. We can thus employ disperionless media 
as perfectly consistent relativistic models to understand 
the key features of some exotic yet realistic effects of light 
in moving dielectrics. 

The optics of moving media has a long history. In 1818 
Fresnel g discovered theoretically that the speed of light 
w in a uniform yet moving medium of refraction index n 
depends on the medium velocity u as 



c 

V = — 






(1) 



So the effective refraction index c/v is changed when the 
medium is moving. In 1851 Fizeau M verified experi- 
mentally Fresnel's motional effect (||). In 1895 Lorentz 
Q derived an additional drag effect that is due to optical 
dispersion. Zeeman [ Ml w as able to measure Lorentz's ef- 
fect. In 1913 Sagnac |[10| observed phase shifts of light in 
a rotating interferometer. In 1925 Michelson, Gale and 
Pearson [|ll|] measured the Sagnac effect of Earth's mo- 
tion in an incredible interferometer by 1925 standards. 
And of course, today's fiber gyroscopes prove that the 
interference of light is sensitive to motion. 

Despite the long history of optics in moving media, a 
sufficiently general theory has been still missing, with 



the exception of Berry's and Klein's parallel work ||] 
and of two earlier papers by Pham Man Qua n [p| , to 
the best of our knowledge. In 1908 Minkowski ]12| pio- 
neered the modern theory of dielectrics. In 1923 Gordon 
published a far-sighted paper M on electromagnetism 
in dielectrics and in gravitational fields. Here he dis- 
covered a deep analogy between gravity and dielectric 
media. According to Antoci and Mihich |13|, Gordon 
1^ also settled the debate about Minkowski's |1^] versus 
Abraham's [[l4| energy-momentum tensor in Abraham's 
favor. However, Gordon considered only very briefly the 
geometrical optics of moving media. Pham Mau Quan 0] 
studied ray propagation in more detail but still not ex- 
haustively. Landau and Lifshitz |]15[ and Van Bladcl p^ 
summarize to some extend the current state of the the- 
ory, but do not focus on the propagation of light. More 
importantly, the motion of the medium has been usually 
assumed to be uniform. Exceptions are the papers p7| ] 
that, however, treat only special cases such as moving 
dielectric boundaries and the consequent modification of 
Snell's law and the Brewster angle. Landau and Lifshitz 
|p5| write explicitly that they "neglect slight effects due 
to the possibility of a velocity gradient" . 

However, these neglected effects are indeed measur- 
able with modern interferometry (see Sec. IV A for an 
estimation). Furthermore, effects due to velocity gradi- 
ents establish interesting connections between the optics 
of moving media and other fields of physics. Hannay 
p8[ discovered an analogy between light in moving me- 
dia and charged matter waves in electromagnetic fields. 
The flow u turned out to play the role of the electro- 
magnetic vector potential. Hannay used path integrals 
in paraxial approximation to arrive at this conclusion. 
Cook, Fearn, and Milonni [|9| analyzed further the con- 
nection between light in moving media and charged mat- 
ter waves, assuming relatively slow medium velocities and 
light that is perpendicularly polarized with respect to 
the flow. The magnetic analogy of light in moving me- 
dia is particularly interesting, because the light propa- 
gation at a fluid vortex corresponds to the Aharonov- 
Bohm effect [^ of electron waves that enclose a local- 
ized magnetic flux. Light that travels through a dielectric 
vortex attains an Aharonov-Bohm phase shift. On the 
other hand, atoms that pass an electromagnetic vortex 
experience an Aharonov-Bohm effect as well [^. In- 
teresting quantization effects arise when the atoms form 
a macroscopic condensate [Q. The magnetic model of 
waves in moving media is not restricted to light. Indeed, 
Berry et al. 123] report both the theory and an experi- 
ment that demonstrates an Aharonov-Bohm effect with 
water waves. Acoustical analogs of the effect have been 
observed in moving classical media [ p4| and are predicted 
for superfluids p5[. However [p|, the magnetic model of 
light in moving media is only valid as long as the medium 
velocities are sufficiently small. In general, the moving 
medium turns out to act rather as a curved metric, i.e. 
as a gravitational field, on the light. Note that Unruh 
lEa arrived at a similar model for non-relativistic sound 



in moving fiuids that also holds for superfluids |^ . 

In Section II we summarize the theory of effects in first 
order of u/c before we turn to the general case in Sec. HI. 
Section II sets the scene by presenting a short review of 
partly published yet not widely known results, whereas 
Sec. HI is the core of our paper. Here we establish the ef- 
fective wave equation, a Hamiltonian, a Lagrangian and 
the metric of light in moving dispersionless media. Ad- 
ditionally, we show how the magnetic model of Sec. II 
and Refs. [^,|9| is related to the gravitational one P,p[ . 
Both Sec. II and Sec. HI derive Fresnel's formula (0)^ at 
least to the lowest order in u/c, seen, however, at each 
case in the light of a distinct physical model. Quantita- 
tive differences between the two concepts are only visible 
in higher order. One example of gravitation-like effects 
is the light deflection at a vortex that we analyze in Sec. 
IV and that resembles the deflection of light due to Sun's 
gravity. In an extreme case the vortex might even appear 
as an optical black hole |^ similar to Unruh's dumb hole 
P6|,p7|. 



II. SLOWLY MOVING MEDIA 

Consider a moving nondispersive dielectric medium 
with refractive index n and flow u. We allow both n 
and u to vary in space and time. However, n and u shall 
not change signiflcantly over the spatial scale of an opti- 
cal wave length and over one optical cycle, respectively. 
In this section we model light waves by a scalar complex 
function -0, for simplicity. In particular, we do not con- 
sider the polarization of light. However, we show in Sec. 
Ill that the propagation of light is indeed independent 
of the polarization, as long as the medium varies only 
gradually compared to optical oscillations. Furthermore 
we assume that the medium moves at moderate velocities 
such that we can restrict ourselves to effects that occur 
within the lowest order in u/c. 

Our starting point is a simple model: Imagine that the 
moving medium consists of small cells or drops. Each 
cell shall be small enough such that the refractive index 
n and the velocity proflle u of the medium does not vary 
signiflcantly. On the other hand, each cell shall be large 
compared to the wave length of light. We thus assume 
that in each cell (in each comoving frame of the medium 
denoted by primes) the optical fleld -0 obeys the wave 
equation 



^"-"iw^y* 







(2) 



An observer sees the light in the laboratory frame. To 
transform the wave equation to the observer's frame, we 
write Eq. (H) as 



V'^- — 






n^- 1 



92 \ 



(3) 



We note that the d'Alembert operator V'^ — d"^ /{cdt'Y 
is a Lorentz invariant, and thus we transform solely the 
remaining time derivatives in the wave equation (^). In 
the lowest order in u/c a temporal change d/dt' in the 
medium frame appears in the laboratory frame as the 
time derivative d/dt plus the local flow u • V. Therefore, 
we obtain in first order 



y'~^ 






- u • V ^ ) V = . (4) 



Note that this derivation of the wave equation for light in 
slowly moving media follows Fresnel's original idea [|6|j28| 
who divided the ether into an invariant part and a second 
part that the medium is able to drag. We will see shortly 
that Fresnel's formula @j is a direct consequence of the 
wave equation (0). 

In the limit of Hamilton's geometrical optics we rep- 
resent the optical field ijj in terms of a slowly varying 
amplitude and a rapidly changing phase. 



-0 = Ae 



iS 



with 



S ~ / (k • dx — Ludt) 



(5) 



(6) 



The wave vector k corresponds to the momentum of a fic- 
titious particle that follows a light ray and the frequency 
u! plays the role of the Hamiltonian. We substitute the 
ansatz (0) and (0) into the wave equation (0) , neglect the 
variation of the amplitude A, and obtain in first order in 
u/c the dispersion relation 



0: 



-OJ 



2uj- 



1 



-u k 






n^ -1 



-uk 



(7) 
(8) 



The Hamiltonian of light rays, H, is equal to the fre- 
quency Lo. We read immediately from Eq. (0) that 



H ^ -k 

n 



1- J_ lu-k 



(9) 



The ray trajectories are solutions of Hamilton's equations 

dx _dH dk _ dH 
dt 9k ' dt dx 



(10) 



The first part of the Hamiltonian, ck/n, describes light 
rays in a medium at rest p9| . The rays avoid regions of 
high refractive index in order to minimize their dimen- 
sionless optical path lengths J k • dx = {i^/c) J n et • dx 
with Gk = k/fc. The second part of the Hamiltonian 
describes Fresnel's drag effect. Indeed, we obtain from 
Hamilton's equations ( |lO| ) 



dx 



-ek 



ek 



k 



(11) 



This is nothing but the vectorial version of Fresnel's orig- 
inal formula (|l]). 

As has been pointed out earlier pt^ , a uniform medium 
in stationary motion acts on light in the same way as a 
magnetic field acts on charged matter waves. To under- 
stand this remarkable connection within our theory of 
ray propagation, we introduce a rescaled ray vector w, 
or, equivalently, a reparameterization of the ray trajec- 
tory, 



w = fcv = — k 

n 



1 T 1 fcu . 



Let us derive an equation of motion for w. 
replace fcu by n{uj / c)\i in first order, and get 



(12) 



First we 



w = fcv = — 
n 



- 1 



-LOVi 



(13) 



Then we use Hamilton's equations ( pO| ) and the relation 

a; du cj , „, / „^ , ^ 

n- — =:?i-(v V)u= (v^f • V)u (14) 

that is valid in first order. We obtain the Lorcntz-type 
equation of motion 



dw 
~dt 



1 



1 



(V X u) X w . 



(15) 



Light rays in slowly moving media behave like charged 
particles in magnetic fields where the flow u ap pears as a 
vector potential. The Lorentz-type force ( [l5| ) conserves 
the modulus of the modified velocity w which is equal to 
Lj in regions where the medium is at rest, 

V? =oj'^ . (16) 

We substitute w by the right-hand side of Eq. (O) and 
retranslate the resulting dispersion relation into a wave 
equation, replacing kf/! by — iVV'- In this way we ob- 
tain exactly the Schrodinger equation of a charged matter 
wave in a magnetic field pfl] 



-iV 



n^- 1 



ip — n —^^ 



(17) 



All these arguments support a magnetic model of light 
propagation in moving media p9[ . The fiow u acts as a 
vector potential that modifies the relation between the 
canonical and the kinetic momentum 



n^-1 



-wu 



(18) 



For example, a rotating rigid glass cylinder will act like a 
homogeneous magnetic field on light that travels inside. 
The rotating cylinder will bend light rays, irrespective of 
their distance from the rotation axis. Another example 
is a vortex that will act like a thin solenoid, see Sec. IV. 
Light rays arc not bent, but, similar to the Aharonov- 
Bohm effect pfl], rays that enclose the vortex attain a 
phase difference M|. 



III. LIGHT IN MOVING MEDIA 

A. Wave optics 

Let us develop a completely relativistic theory of light 
propagation in moving nondispersive media. Like in Sec. 
II we assume that the refractive index n and the flow u 
do not vary significantly over one optical wave length and 
one optical cycle, respectively. We neglect the dispersion 
of light, i.e. the frequency dependence of n. We allow ar- 
bitrary medium velocities and we will employ a covariant 
notation Q. Our starting point is the following postu- 
late: In all locally comoving medium frames (denoted 
by primes) the electromagnetic field-strength tensor F' 
shall obey the wave equation 



V'- 



.2 Qf2j M- 







(19) 



Note that this postulate uses implicitly the assumption 
that the refractive index varies only gradually. Otherwise 
additional terms become important in the wave equation 
| p5[ , terms that describe polarization changes (at the sur- 
faces of dielectrics, for instance). 

Let us transform the wave equation (03) to the labo- 
ratory frame. As a first step we reformulate Eq. (|lj) in 
a covariant notation. We employ the four-gradients 



a.' 



cdf 



v 

,v 






i_a_ 

cdt' 
ld_ 



-V 

-V 



and the four-vector field of the medium flow 



w = 7 



V c 



Ui/ = 7 



V c 



with the relativistic factor 



7- 1 



2\ -1/2 



(20) 



(21) 



(22) 



In a comoving medium frame the four-vector u"^ is lo- 
cally 



(23) 



^"^ = (1,0) 



Therefore, we can easily write our starting point ( |19| ) in 
the covariant expression 



[d'^d"^ + in^-l){u"^d'J']F'^^^0 



(24) 



Throughout this paper we employ Einstein's summation 
convention. When we transform the wave equation (EJ) 
to the laboratory frame we should transform both the 
derivatives and velocities, and the field-strength ten- 
sor. A Lorentz transformation of a tensor depends of 
course on the velocity of the moving frame [pl[. Since 



the medium velocities vary only gradually compared with 
the rapid oscillations of F', we can neglect the deriva- 
tives of the Lorentz transformations of F' in Eq. (|2J). 
In other words, the wave equation ( |24|) is valid both 
for F' and i^^,y. Furthermore, the differential operator 
d'^d'°' + {n^ — l){u"^d'^)'^ is a Lorentz scalar, and therefore 
we can simply drop the primes in Eq. (|24|), to arrive at 



[a„9" + (n2-l)(M"a„)2]F^, = 



(25) 



This wave equation describes the propagation of light in 
a moving nondispersive medium, provided that both the 
refractive index n and the fiow u varies only gradually. 

We see that the final wave equation ( p5| ) does not mix 
the components of the field-strength tensor. The prop- 
agation of light beams does not depend on the polariza- 
tion, i.e. moving media are not birefringent. This result 
justifies the scalar model of Sec. II. Note however that 
the transport of the field amplitudes along light beams 
is certainly polarization-dependent. To describe this ef- 
fect one should consider, instead of the wave equation 
(PSI ) , the complete set of Maxwell's equations in moving 
dielectrics 0Jl5|,ra. 



B. Geometrical optics (Hamiltonian) 

How does a moving medium act on light rays? Accord- 
ing to Hamilton's geometrical optics, we try the eikonal 
ansatz iBll 



^ nu — ^ tiu G ~r C.C 



iS 



(26) 



with 



S" = f [k-dx- udt) = - I Kdx" . (27) 

Here we have employed the four-differential 

dx" ^{cdt,dTi) (28) 

and the wave four-vector 

fc, = (-,-k) =-9,5. (29) 

Assuming a rapidly changing phase S compared with the 
envelope .F^jy, we derive from the wave equation (Gq) the 
Hamilton-Jacobi equation of light rays. 



with 

using the flat Minkowski metric 



V''" = V^. 



1 
1 



(30) 
(31) 

(32) 



Explicitly, we get the dispersion relation 



(33) 



To find a Hamiltonian for light beams we solve Eq. 
(p3|) ior Lu = H. We obtain 



H = 



1/2 



2; 2 

C k 



r(u-k)' 



rUk 



1/2 



(34) 



In first order in u/c we recognize our previous result (S). 
Ray trajectories follow from the Hamiltonian (|4|) as so- 
lutions of Hamilton's equations (p^. Here the time t 
plays merely the role of a parameter to describe the tra- 
jectories. Of course, infinitely many parameterizations 
exist that result in equivalent ray trajectories but stem 
from different Hamiltonians. Does our particular param- 
eterization have a physical meaning? Let us consider the 
formal velocity v' in the comoving medium frame llq] , 



V — u + (7 — l)u(u • V — u^)/u^ 
7 (1 — u • v/c^) 



(35) 



We use the first of Hamilton's equations ( jlO| ) and the 
Hamiltonian (p^), and obtain after some algebra 



(36) 



Light travels with the velocity of light. Therefore, we 
can identify the formal ray parameter t with the physical 
travel time of light in the laboratory frame. 

Apart from elucidating the physical meaning of time 
for our Hamiltonian, the relation (|3^) is simply an ex- 
plicit conservation law during the ray propagation. Con- 
servation laws are connected to symmetries. In fact, the 
Hamiltonian (B4) has the remarkable structure 



H = k h{() , C = u • ek , ek = - 



Consequently, 
_ dH 



= /iGk 



dh 



Gk X (U X Gk) 



(37) 



(38) 



We see that the velocity vector v is independent of the 
wave number k. The trajectories of light beams in mov- 
ing media do not depend on the wave properties of light. 
Geometrical optics involves indeed solely the geometrical 
aspects of light rays. 



C. Lagrangian 

We have studied the propagation of light in moving me- 
dia in the spirit of Hamilton's geometrical optics. Can 



we find a Lagrangian as well? Let us calculate the La- 
grangian directly from the Hamiltonian, using the struc- 
ture (^^ and (|8|) of the Hamiltonian and the velocity 
vector, respectively. We get 



--f- 







(39) 



The Lagrangian vanishes. Furthermore, we cannot ex- 
press the canonical momentum k in terms of the velocity 
V, because, according to Eq. (p8|), v does not depend on 
the modulus of k. The encountered problems in intro- 
ducing a Lagrangian for light rays in moving media have 
been known for the special case of light in vacuum |3l| ] 
and have been attributed pU to the vanishing rest mass 
of light. 

Let us provide light with an artificial rest mass m that 
we let approach zero at a later, appropriate moment. We 
replace the right-hand side of the wave equation (p5|) by 



m c F^^^ and we obtain in the limit of geometrical optics 
the Hamilton-Jacobi equation 



g'"'{d^S){d,S)^m'^ 



(40) 



We wish to find a Lagrangian representation of the phase 
S such that 



Ldt 



For this we introduce the matrix 



1 



Qtiy = 7?Mi' + ^ ^ 1 "m" 



(41) 



(42) 



One verifies easily that g^v is the inverse of of g'*'', uti- 
lizing that Uvu'^ is normalized to unity. 



5^05"" = 5^^. 



(43) 



We obtain the solution of the Hamilton-Jacobi equation 

S = -mc I ds (44) 

with 

ds^ = g^.^dx^'dx" . (45) 

To prove this result we note that the wave four-vector 
gives 

dx"^ 
kfj. = -d^S = mcg^„—j— , (46) 

and thus satisfies the Hamilton-Jacobi equation (Eoh 



ds ds 



(47) 



The expression (M) of the phase S has the desired La- 
grangian structure. 



dxt^ dx'^ 



'^'-^nr'^' 



S — —rac 
with the exphcit Lagrangian 



Ldt 



L 



-mc \ c-^ — v^ + 



1 



1 7^ c 



(48) 



(49) 



D. Metric 

According to the action principle, hght rays minimize 
the phase S. Because the phase (Q) is proportional to 
the artificial mass m, the minima of S do not depend 
on the numerical value of m at all. We thus arrive at 
an entirely geometrical picture of ray optics in moving 
media, in complete analogy to the equivalence principle 
of general relativity. Light rays are geodesic lines with 
respect to Gordon's metric M 



ds^ = T^f.^dx^dx" + ( ^ - 1 ) [u^dx'^y 



V^u^ 



1 c^dt'^ 



rdt" - dx'- 



(50) 



(51) 



introduce the contravariant wave vector with respect to 
the metric of the moving medium, 



fc^ EE 5^"/^^ 



(55) 



According to Eqs. (^) and ( p^ ) the fc^ four-vector is 
proportional to the four-velocity. 



k^' ^ 



dx^' 



dt dxf' 



(56) 



ds ds dt 

Because dx^/dt is equal to c, by definition (|2l), we obtain 



k" 



k° dxt" 
c dt 



(57) 



a relation that remains valid for vanishing m and ds. In 
the magnetic model of light propagation in moving me- 
dia, see Sec. II, the three-dimensional k vector differs 
from the velocity by an additional term that is propor- 
tional to the flow and that plays the role of a vector 
potential. However, as we have seen, the contravariant 
fc^ vector is in fact proportional to the velocity dx'^/dt. 
Therefore, the appearance of the flow as a vector poten- 
tial is an illusion that we get when we do not discriminate 
between co- and contravariant vectors, without yet ap- 
preciating the inherent space-time geometry of geomet- 
rical optics of moving media. 



While traveling in a moving medium, light minimizes the 
proper time measured with respect to c/n. And, because 
light actually travels with the velocity of light ( Pq ) , the 
proper time is zero on the ray trajectory. Light rays 
follow zero-geodesic lines of the metric (pO|). 

A moving medium acts as a curved space-time on light. 
In terms of Riemann's geometry, a medium with uniform 
refractive index n turns out to have the scalar curvature 

M 



R 



{n^ - If 
2n2 



{u^'u'' -l^^'''){^^,u'){^,u,) 



+ (1 - n^) {2u''d,d,,u^ + {d^un^) 



(52) 



as one obtains after some algebra. When the flow obeys 
the continuity equation. 



d^u"" = , 
the curvature simplifies to 



(53) 



R = 



2vP- 



(w^u''-77'^'')(a^«')(a.M,) 



+ ^;^(5X)(5.«'')- 



(54) 



In curved space-time the distinction between co- and 
contravariant objects is particularly profound. Let us 



E. Slowly moving media 

In Sec. II we summarized the magnetic model [|l8[|l9[ 
of light in slowly moving media. Let us formulate the 
first-order effects in terms of the geometrical model [|] . 
To the lowest order in u/c we obtain the metric 



ds" 



[cdtf - dx^ + 2 ■ 



-1 u 



• dx cdt . 



(58) 



Therefore, the covariant metrical tensor g^j^n is given by 
the matrix 



1 



- 1 u 



9^^u 



^2-l u 



(59) 



To find the contravariant metrical tensor we compare the 
dispersion relation (0) with Eqs. (|29| ) and (|30|), and read 
off the result 



^M^ 



in' 



(r? 




(60) 



In the limit of low medium velocities, three-dimensional 
space appears to be flat to light that travels in the 
medium, yet four-dimensional space-time is curved. Let 
us compare co- and contravariant wave vectors. As a 
consequence of Eqs. (|l^) and (|9|) and of the definition 



( p5| ) , the space components fc' of the contravariant fc^ are 
given by 



k+{n^ 



LO n 

1) — u = — w 

c c 



— fcv 
c 



(61) 



Therefore, the contravariant wave vector appears as the 
kinetic momentum, in contrast to the canonical one that 
is represented by the covariant wave vector. Finally, 
we solve Eq. (p^) for k/fc and arrive at the vectorial 
form (^l]) of Fresnel's classic result (||), seen here in the 
Riemann-geometrical model of light in moving media. 



IV. LIGHT AROUND A VORTEX 

A. Optical Aharonov— Bohm effect 

Consider a vortex flow in an incompressible liquid. The 
vortex may be created by the action of a mixer or simply 
by letting the liquid flow out off the container. To a good 
approximation the velocity profile of the vortex is given, 
using cylindrical coordinates, by the expression p3] 



u = — e,, 



(62) 



Let light travel through the whirling liquid. In first order 
of m/c, light experiences the medium flow in the same way 
as a charged matter wave experiences a vector potential 
|Qq,19[. In particular, the vortex ( p2| ) corresponds to an 
infinitely thin solenoid that generates a strong magnetic 
field inside yet no field outside. 



V X u = 



(63) 



However, as Aharonov and Bohm discovered in their sem- 
inal 1959 paper |2^, a charged matter wave will attain 
a phase shift without feeling a force, and so will light. 
We compare the Schrodinger equation of Aharonov and 
Bohm 1221 with the wave equation (|l^ and read off the 
optical Aharonov-Bohm phase shift 



'Pa 



— 2-KV. 
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(64) 



How large is the effect? For estimating the order of mag- 
nitude we assume a refractive index n of 1.5 and an opti- 
cal frequency lo of 3 x lO^^s"^. Furthermore we assume 
that the liquid circulates at a radius r of 2 cm with a 
frequency u/r of 10 s~^. We obtain a phase shift (/j^^ 
of 10~'^. The effect is small, yet one could enhance it 
significantly by letting the light travel through the liq- 
uid many times in an interferometer. In this case the 
phase shift (/^^^ is effectively multiplied by the number 
of round-trips. One could also perceive the interferome- 
ter as a resonator where the effect of the moving medium 
changes the spatial mode profile. The precision of mod- 
ern interferometry is certainly sufficient to detect the o 
tical Aharonov-Bohm effect. (Remember that Fizeau 



has seen the precursor of the effect as early as 1851.) Fig- 
ure 1 illustrates the long-range behavior of light waves 
that pass a vortex flow. 
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FIG. 1. Light wave passing through a vortex flow. Picture 
(A) illustrates the optical Aharonov-Bohm effect. Light is 
incident from the right, and we plot the real part of the op- 
tical fleld Tp that is subject to the wave equation (Uj). We 
use for ip the well-known Bessel series due to Aharonov and 
Bohm [20]. Picture (B) sketches the wave patterns of light, 
taken into account the slight ray bending that we discuss in 
Sec. IV C. We use the same scale as in (A). Before reach- 
ing the vortex the wave fronts bend and then they split at 
the vortex core. Behind the core the two parts of the wave 
interfere. We created the picture by following a bundle of 
exact ray trajectories IjPfX)- A wave front is defined by the 
points reached at a particular value of the time parameter. 
Although this approach does not provide information about 
the wave amplitudes, the picture is quite accurate, because 
mostly only one part of the split wave contributes to interfer- 
ence, except in regions very close to the cut at the negative 
X axis. We see that picture (B) reproduces the prominent 
long-range features of the Aharonov-Bohm wave depicted in 
(A). 



the relativistic Euler equation. For this we construct the 
energy-momentum tensor [ p3| 



and assume the pressure profile 

>V2 



p = pc exp 
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(69) 



(70) 



Here p denotes the constant mass density of the incom- 
pressible fluid. We obtain by direct calculation that T'"' 
satisfies the Euler equation 



d^T^"" = 



(71) 



This proves that the vortex ( pq ) is indeed a possible flow 
that, as we have seen as well, generates the pressure pro- 
file (ItO). At sufficiently large distances from the vortex 
core the pressure behaves like 



B. Relativistic vortex 

Let us analyze vortex effects to higher order in u/c. 
Note however that the non-relativistic vortex (p2) per- 
mits medium velocities that exceed the speed of light, 
when taken seriously near the vortex core where higher- 
order effects are to be expected. Let us therefore seek a 
proper relativistic vortex. We anticipate that the flow is 
still circular. 
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(72) 



The rapidly falling pressure will attract the surface of 
the liquid and create a hole at the vortex core. Let us 
nevertheless treat the hole as being infinitely thin. Near 
the vortex core the pressure is finite 



lim p — —pc 



(73) 



and the medium velocity reaches the speed of light 



u{r)e^ 



(65) 



but that we must correct the non-relativistic velocity 
profile (|6|) by the relativistic 7 factor (gg), 






(66) 



to prevent the medium from becoming tachyonic. Given 
the velocity profile (|66|), we solve Eq. ( p^ ) for 7, and 
obtain 



lim u = c 
r— >o 



(74) 



At large distances from the core the medium flows as 
described by the non-relativistic Eq. (p2). Of course, 
dramatically relativistic vortices are not to be expected 
to be created in experiments on Earth in the near future, 
but they might exist as astronomical objects. 



C. Light bending 



7= Wl 



W2 



Therefore, the four-flow (^) is finally 
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Let us study the propagation of light rays around the 
relativistic vortex ( p a t) - For this we solve the Hamilton- 
Jacobi equation ([lOPthat reads explicitly 



s-" - c^ivsy + (n^ - i)Y (s + uvsy 







(75) 



We obtain for any circular flow (pa) the solution 



Of course, so far we have simply guessed the proper rel- 
ativistic generalization of the vortex flow (p2h. To prove 
that our result is indeed physically meaningful, we shall 
show that the four-flow (pq) is consistent with relativistic 
hydrodynamics [ p3[ . First we note that any circular flow 
( p5| ) satisfles the equation of continuation ( p3| ) . Then we 
prove that the particular four-flow (pq) is a solution of 
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Let us first study light rays that travel sufficiently far 
outside the vortex core. We use our velocity profile (|6q ) 
with the factor (67), expand {dR/dr)^ to quadratic order 
in r~^, and get 



dR 

1^ 



with 



I, 






l + in^-1)^ 
c 



n'W 



(78) 



(79) 



As we shall see below, Z^^ plays the role of the kinetic an- 
gular momentum that is altered by the optical Aharonov- 
Bohm effect. According to canonical mechanics |Q tra- 
jectories are given requiring that the derivative of the 
action with respect to generalized canonical momenta is 
a set of constants. In particular we put 



dS_ 
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In this way we obtain the trajectories 
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with 
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and, after solving the integral, 
I 
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To first order in u/c we disregard the W^/c^ term in Eq. 
(H) and get from Eq. ^ 



rcos{ip - ipo) = ro ■ 



(84) 
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to second order in W/c. Similar to the bending of star 
light due to Sun's gravity |^, hght is attracted by a 
vortex and is consequently slightly deflected. We express 
ifj^ in terms of the ffi'st-order phase shift (p5), and obtain 
to second order in W/c 
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The deflection angle is extremely small for small 
Aharonov-Bohm phases 27rz/^g . Figure 2 illustrates the 
bending of light at a vortex. The figure shows also light 
rays in the immediate vicinity of the core, a case that we 
shall consider in the next subsection. 




Therefore, far outside the vortex core light travels along 
straight lines (in accordance with our previous consid- 
erations on the optical Aharonov-Bohm effect). Let us 
study indications of higher-order effects. At infinity the 
angle (p approaches 

TT I 

(^oo = lim ip^ipo + -^^ . (85) 

r^oo z nrQ 

The impact parameter of the ray is given by 

b = lim rsm{ip - >f^) ^ — — . (86) 

r — >oo 77, 

This result shows explicitly that l^^ corresponds indeed 
to the kinetic angular momentum. When the light leaves 
the vortex region the ray is deflected by the angle |34 




FIG. 2. Trajectories of light rays in the vicinity of a vortex 
flow. Each picture shows a pair of trajectories belonging to 
the same value of the angular momentum. In one case the 
light ray approaches the vortex from infinity, is bent, and dis- 
appears to infinity. In the other case the light stems from the 
origin and falls back after having reached its maximal range 
from the vortex core. At distances not covered by these tra- 
jectories no light rays are allowed that belong to the chosen 
particular value of the angular momentum. (A) Negative an- 
gular momentum (positive impact parameter). Light swims 
against the current. (B) Positive angular momentum (nega- 
tive impact parameter). Light travels with the flow. 



D. Optical black hole 

A vortex attracts light like any other test particle, be- 
cause of the rapidly falling pressure profile (JT^). Can 
light fall into the vortex core? To answer this question 
we study the turning points of light rays, using the com- 
plete relativistic action rtTq) with the radial component 
([TTI). Turning points are the zeros of (dR/dr)^ because 
beyond these points the action S is purely imaginary and 
hence covers a forbidden region. We represent Eq. (177 
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We see that each radius r can be reached as a turning 
point of two trajectories that are characterized by the 
angular momenta 



nc:Lu 
nu ± c 



(90) 



The point may be an inner or an outer turning point, 
depending on the sign of 



c± 



d fdR 



dr V dr 



(91) 



The vortex may emit light that begins to fall back to 
the core at an inner turning point with negative C. Or, 
alternatively, incident light from outside the vortex is 
bent and comes to the closest distance to the core at 
an outer turning point with positive C,. To distinguish 
inner and outer turning points we calculate C± and write 
the resulting expression in terms of the derivative of l± 
regarded as a function of r. We obtain 



dr ^ dr r (nu ± c)^ 



(92) 



The velocity profile u{r) of the vortex (pq) is monotoni- 
cally decreasing and positive. (We assume a positive W. 
Otherwise one could invert the system of coordinates to 
arrive at a positive vorticity.) Therefore, 1+ is positive 



and also yields a positive C-i-- Given a positive angular 
momentum, any point can be reached as an outer turning 
point of incident light rays. Let us turn to the /_. The 
angular momentum L of a turning point is negative out- 
side the radius where the medium velocity reaches the 
speed of light in the medium, c/n. Inside this radius 
the trajectory that corresponds to L approaches an in- 
ner turning point, because C.- is negative. On the other 
hand, for negative angular momenta a critical radius r^ 
exists where (- vanishes. Beyond this radius all light rays 
with negative I must fall into the vortex core. The vortex 
appears as a black hole to light with negative angular mo- 
menta. Let us calculate the optical Schwarzschild radius 
r,. We utilize that 



ru' ~ —u 1 



(93) 



for the relativistic vortex (B8|) and arrive at a third-order 
equation for the velocity Us at the Schwarzschild radius 
where dl-/dr vanishes. 
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This equation has three real solutions labeled by i G 
{1,2,3}, 
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(95) 



The values of the physical solution must lie in the inter- 
val < Wi < c in order to be velocities allowed by the 
theory of relativity. The refractive index n is confined 
to the interval 1 < n < oo. The only function assuming 
physically permitted values at the ends of this interval is 
U'iin) with 



«3(1) 



lim u-i(ri) = 



(96) 



In order to prove that the results for intermediate values 
of n are physical as well, we show that U3 (n) is a mono- 
tonically decreasing function. The first term in (|9g) is 
a product of two positive and monotonically decreasing 
functions and thus is a positive monotonically decreasing 
function itselL So is the second term that is subtracted. 
In order to show that the total expression stays positive 
one can employ the estimates 



1 



arccos(x) > — Vl + xn + 



for 



1< a; <0 



cos(x) > —X 



for — < X < TT. 

2 - - 



(97) 



This leads to an algebraic inequality that can be solved 
in a standard way, thus proving our assumption. In an 
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analogous way one can show that mi and M2 do not as- 
sume physical values. Finally, to calculate the optical 
Schwarzschild radius r^., we solve Eqs. (pw and (pTh for 
r, and get 



B 



w 
U3 



(98) 



Optical black holes could be created using highly dis- 
persive dielectrics [Q,^. These media are distinguished 
by an extremely low group velocity of light. To some 
extend, they resemble dispersionless dielectrics with an 
extremely large refractive index. Let us therefore con- 
sider the limit of large n. We obtain from Eq. (^5|) the 
critical flow velocity 



"3 



c 

2n 



(99) 



A light ray with a negative angular momentum, i.e. light 
that swims against the current, is trapped when the 
flow velocity approaches half of the speed of light in the 
medium. Figure 3 illustrates the fall into the optical 
black hole. 





FIG. 3. Optical black hole. The pictures show trajectories 
of light rays falling into the vortex core. In case the angular 
momentum of an incident ray lies within a critical interval, 
the ray is dragged by the vortex flow and flnally falls into the 
singularity. (A) We have chosen several values of the angular 
momentum that cover most of the critical interval. (B) The 
chosen values of the angular momentum lie immediately above 
the negative bound of the critical interval. In this case a 
trajectory is very sensitive to small changes in the angular 
momentum. The radius of the limiting circle is the optical 
Schwarzschild radius. 



E. Scaling 

A vortex flow can cause an optical Aharonov-Bohm 
effect and may force light to fall into an optical black 
hole. Note that these phenomena depend, in principle, 
entirely on the mere presence of the vortex and not on 
the particular value of the vorticity, even when W is very 
small. To see this we study the scaling properties of the 
relativistic vortex. A closer inspection to Eqs. (p5m^) re- 
veals that the dimensionless velocity profile u/c depends 
only on the combination (c/W) r. Let us introduce the 
new variables 



t=—t,r^—r,l^—l,S^—S. (100) 

Due to the scaling of u/c the scaled eikonal S* is a so- 
lution of t he Hamilton-Jacobi equation (^) in the vari- 
ables (100). The scale of I and S was chosen such that S 
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corresponds to the explicit solution ([T^) with the param- 
eters uj and I. Note that the Hamilton- Jacobi equation 
( [75| ) in the dimensionless coordinates t, f, and (p does 
not contain W anymore. Consequently, the relativistic 
effects on light at a vortex flow do not disappear when 
the vorticity approaches zero. However, the smaller the 
vorticity is the higher should be the frequency of light, 
in order to produce an equivalent Aharonov-Bohm phase 
shift (|6^). This is caused by the scaling of the eikonal S 
that governs the phase profile of the optical field. Fur- 
thermore, for a low vorticity the optical Schwarzschild 
radius is accordingly small. In our idealized model, any 
vortex is an optical black hole, irrespective of the value 
of W. In practice, the incident light is more likely to hit 
an obstacle that we have entirely ignored — the vortex 
core. Similar to a star that acts as a black hole when 
the gravitational Schwarzschild radius exceeds the star's 
size, a vortex is an optical black hole only if the core is 
smaller than the optical Schwarzschild radius. 



V. SUMMARY 
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A moving dielectric appears to light as an effective 
gravitational field [|| . At low fiow velocities the dielectric 
acts on light in the same way as a magnetic field acts on 
a charged matter wave [18 19 1. The fiow plays the role 
of the vector potential. We have shown how the two ef- 
fective models are related to each other: The covariant 
wave vector corresponds to the canonical momentum of 
a light ray, whereas the contravariant wave vector plays 
the role of the kinetic momentum. In curved space-time, 
and hence in a dielectric, co- and contravariant vectors 
are distinct. For low medium velocities they differ in 
precisely the same way as the canonical momentum of a 
charged particle differs from the kinetic one by the vec- 
tor potential. Additionally, we have derived and studied 
an effective Hamiltonian for ray propagation in moving 
media. Within the limits of geometrical optics, ray tra- 
jectories do not depend on the polarization. Introduc- 
ing a fictitious rest mass for light we have constructed a 
corresponding Lagrangian. This way has lead us to the 
description of a medium in terms of a metric [^1 . 

First-order optical phenomena are certainly detectable 
using modern interferometry and ordinary dielectric flu- 
ids. In this way one could infer the velocity profile of 
transparent incompressible liquids from interferometric 
measurements. To observe some spectacular relativistic 
effects, one could employ highly dispersive quantum di- 
electrics [Q . One can demonstrate an optical Aharonov- 
Bohm effect and create an optical black hole with a quan- 
tum vortex as the center of attraction ||5[] . Our paper has 
established a consistent yet simplified model of such phe- 
nomena. We have thus reason to hope that we may stim- 
ulate experiments to demonstrate gravitational effects on 
Earth that usually belong to the realm of astronomy. 
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